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I consider the hydrodynamic stability of imploding ideal gases as an idealized model 
for inertial confinement fusion capsules, sonoluminescent bubbles and the gravitational 
collapse of astrophysical gases. For oblate modes (short-wavelength incompressive modes 
elongated in the direction of the mean flow), a second-order ordinary differential equation 
is derived that can be used to assess the stability of any time-dependent flow with planar, 
cylindrical or spherical symmetry. Upon further restricting the analysis to homologous 
flows, it is shown that a monatomic gas is governed by the Schwarzschild criterion for 
buoyant stability. Under buoyantly unstable conditions, both entropy and vorticity fluctu¬ 
ations experience power-law growth in time, with a growth rate that depends upon mean 
flow gradients and, in the absence of dissipative effects, is independent of mode num¬ 
ber. If the flow accelerates throughout the implosion, oblate modes amplify by a factor 
(2C')l Ar °l ti , where C is the convergence ratio of the implosion, No is the initial buoyancy 
frequency and ti is the implosion time scale. If, instead, the implosion consists of a coast¬ 
ing phase followed by stagnation, oblate modes amplify by a factor exp (tt | A r o| t s ), where 
Nq is the buoyancy frequency at stagnation and t s is the stagnation time scale. Even 
under stable conditions, vorticity fluctuations grow due to the conservation of angular 
momentum as the gas is compressed. For non-monatomic gases, this additional growth 
due to compression results in weak oscillatory growth under conditions that would oth¬ 
erwise be buoyantly stable; this over-stability is consistent with the conservation of wave 
action in the fluid frame. The above analytical results are verified by evolving the com¬ 
plete set of linear equations as an initial value problem, and it is demonstrated that 
oblate modes are the fastest-growing modes and that high mode numbers are required 
to reach this limit (Legendre mode i > 100 for spherical flows). Finally, comparisons are 
made with a Lagrangian hydrodynamics code, and it is found that a numerical resolu¬ 
tion of ^30 zones per wavelength is required to capture these solutions accurately. This 
translates to an angular resolution of or < 0.1° to resolve the fastest-growing 

modes. 
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1. Introduction 


Ex amples of imploding gases i nclude inertial confinement fusion (ICF) capsules (jAtzeni fc Mever-ter-Vehn 
2004h . sonolu minescent bubbles ( Suslick fc Flanniganll2008 ) and core-collapse supernovae 


djankal 2012 1. The hydrodynamic stability of these flows is an important issue, as pertur¬ 
bations of sufficient amplitude can drain the energy driving the implosion, and the break¬ 
down of symmetry, even for small amplitudes, can have important diagnostic/observational 
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effects (e.g. iMurnhvi 12014 1. Stability analyses of an imploding gas are significantly com¬ 
plicated by the time-dependent mean flow, and homologous flow (i.e. flow in which fluid 
elements share a common time-dependent scaling, the Hubble flow being a prime ex¬ 
ample) is a useful idealization that allows some analy tica l progress to be made . Such 
a study was performed for core-collapse supernovae by Goldreic h fe We beil dl980), wh o 
claim stability, but th eir results h ave recently been revisit ed by Cao &; Lou ( 2009 . 201dh . 
who claim instability. Chu ( 1996h applied the analysis of Goldreich fe WebeiT i 1980h to 
a sonoluminescing bubble and also claimed stability. The extensive body of work on the 
stability of ICF implosions has focused almost exclusively on inter faci a l instab ility, with 
little attention being given to the stability of the ga s (see, however. I Greenspan fe Bennev 
1 1 Qfi.'ll : iMiolsness fc Runnell 1978 : Gook et oil l200flh . At the same time, hot-spot turbu¬ 
lence has recently received attention as a potential source of yield degradation in ICF 
capsules ( Thomas fc Kares 2012 : Gatu et al. 20131; Clark et al. 2013 :|Cerian et aZI 20131. 


and t he origin of these vortical flows, if present, remains unclear. [Basko fc Murakami 
( 19981) considered the possibility of buoyancy instability in self-similar implosions but 
did not perform a formal stability analysis. 

The purpose of this work is to clarify some of these issues by performing a sta- 
bilit y an aly sis that co mplements previous work. By solving the initial value problem 
( Lai fc Goldreich 200(j h rather than decomposing perturbations into radial modes, I am 
able to obtain results that are more physically transparent than those of previous au¬ 
thors, including a precise stability criterion. The following assumptions are made: (i) the 
implosion is externally driven rather than driven by gravity; (ii) dissipative effects are 
ignored (this is almost certainly unrealistic for ICF capsules and bubbles, but it permits 
isolation of the instability driving mechanism); and (iii) perturbations are assumed to be 
short-wavelength, incompressive and elongated in the direction of the implosion. I shall 
refer to this final approximation as the oblate limit, and it is perhaps the most useful 
new result to come out of this work, as it allows for significant analytical progress to 
be made in a problem that must otherwise be treated numerically. Under these assump¬ 
tions, I demonstrate that homologously imploding i dea l gases are essen tially governed 
by the Schwarzschild criterion for buoyant stability ( Schwarzschildl 1992 1. with a slight 
modification due to compression. 

The basic equations and mean flow are outlined in gj the stability analysis is given in 
(J31 and g] summarizes and discusses applications. 


2. Basic equations and mean flow 

The fundamental equations used are the continuity equation, Euler’s equation, and the 
first law of thermodynamics for adiabatic flow: 


dp 

dt 


+ pV ■ v = 0 , p 


dv ds 

7t +Vp = ,> ' S = 0 ' 


( 2 . 1 ) 


Here p, v, p and s are the mass density, velocity, pressure and entropy of the gas 
and d/dt = d/dt + v ■ V is the Lagrangian derivative following a fluid clement (e.g. 
[Landau fc Lifshit3 1987 ). I will assume throughout an ideal-gas equation of state, p = 
(7 — 1) pC v T , where 7 , C v and T are the adiabatic index, specific heat and temperature 
of the gas. For this equation of state, s = In (pp -7 ) to within a constant factor; this 
expression will be used in what follows to define the entropy. 

Fluid elements in a homologous flow obey the relationship r(t,a) = r(0,a)h(t) = 
aro\h{t), where r is the spatial coordinate in the direction of the mean flow and the 
Lagrangian label for a fluid element is a = r(0,a)/roi- For a spherical flow, r(0,a) is 
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the field of radial positions for all the fluid elements at t = 0 and roi = r(0, 1) is the 
outer radius of the gas at t = 0; this implies that the scale factor h(t) is normalized 
to unity at t = 0, i.e. h( 0) = 1. The velocity field associated with homologous flow is 
v r (t,a) = aro\h(t) = r(t, a)h/h, where an overdot denotes a time derivative. Mass and 
entropy conservation under spherical adiabatic homologous flow imply p(t, a) = po(a)/h 3 
and p(t,a) = po(a)/h 37 , where a zero subscript denotes a spatial profile at t = 0 , and 
Euler’s equation becomes 


h^~ 2 h=-- 


1 dp 0 


Por-Qic 


da 


= ± 


tv 


( 2 . 2 ) 


where the plus (minus) sign is associated with a decelerating (accelerating) implosion 
and t c is a characteristic time scale. 

For 7 = 5/3 and ho = 0, t he time-dependent portion of expression (12721) yields the 
well-known Kidder implosions ( Kidder 1974 : Atzeni fc Mever-ter-Vehn 20041) 


= \/!± ( J 


(2.3) 


For h 0 ^ 0, an additional solution to (12.21) valid for any 7 is ( Goldreich &; Weber 1980h 

! _ 

-1 

(2.4) 


M0 = 1 - 


'[37-if t 


6 [7 - 1 ] t c 


This solution, while valid only for accelerati ng implosions, is useful for asse ssing the 
impact of the gas equation of state. Following I Atzeni fe Mever-ter-Vehn ( 2004 ). acceler¬ 
ating solutions are defined as starting out from t = 0 , so that they evolve from h = 1 
to h = C -1 , where C = initial Afinai is the convergence factor, i.e, the ratio of initial to 
final sizes. The decelerating solution, on the other hand (the Kidder implosion with the 
+ sign), is defined as stagnating at t = 0, so that its implosion phase evolves from h = C 
to h = 1 . 

Aside from the equation of state, the space-dependent portion of equation ( 12 . 21 ) is 
the only constraint on the spatial profile of the mean gas quantities. This allows for a 
significant amount of flexibility in setting up various mean flow profiles. An analytical 
model that mimics the hot spot of an ICF implosion during deceleration (this is likely to 
be a good approximation for a collapsing bubble as well) is given by 


Po = Po P e exp 1 - e 2 „: 


T 0 = T 0p e , po = pop exp 1 - e 2 


(2.5) 


( 2 . 6 ) 


where the gas quantities have been normalized to their value at a p (the location of the 
density peak; see Fig. [[]). Both a p and a can be regarded as free parameters of the 
model, and can be alternatively expressed in terms of ratios of physical quantities. The 
characteristic time scale in this case is t c = \/ 7 ?'oi cr /cop, where cq p is the sound speed at 
a p . The square of the Brunt-Vaisala, or buoyancy, frequency, 

1 dpds 
7 p dr dr ’ 

is useful for assessing buoyant stability ( Johnson fc Gammi3 2005 ). For homologous flow, 


N 2 = -- 
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Figure 1. Profiles of density (solid), temperature (dotted), pressure (dashed) and the negative 
portion of Nq (dot-dashed) for expressions (12.51) (12.71) with o = 0.2 and a p = 0.9. All quantities 
have been normalized to their peak magnitude. 


N 2 = NqH 1 3t , and for the hot-spot configuration described above, 


2.2 


m 



Mean flow profiles for the hot-spot configuration are shown in Fig. [T] 


(2.7) 


3. Stability analysis 

For an equilibrium configuration (v = 0), it is well known that the local dispersion 
relation for short-wavelength incompressive modes governed by equations ED is 


UJ 


2 



(3.1) 


where ui is the wave frequency and k \ and k r are wa venumbers perpendicular and parallel 
to the mean gradients (see e.g. ljohnson fc Gammia l2005n. It can be seen from expression 
ED that N 2 < 0 denotes ins tability; N 2 > 0 is the Schwarzschild criterion for buoyant 
stability ( Schwarzschildl 1992b . The buoyancy frequency can be rewritten as 


N 2 = ^ds = _ (9 _ g fldT 7 -l d P \ 

7 dr 9 V c s p dr) 7 \T dr p dr) 1 


where g = dv r /dt is the acceleration and c s is the mean sound speed. It can be seen 
from this expression that configurations with “order-over-disorder” (ds/dr < 0) are 
buoyantly unstable. If density gradients dominate, a “heavy-over-light” configuration is 
unstable (dp/dr > 0), whereas if temperature gradients dominate, a “cold-over-hot” 
configuration is unstable (dT/dr < 0). The former situation represents Rayleigh-Taylor 
instability, and the latter represents thermal convection; in general, buoyancy instability 
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depends upon the mean entropy gradient. Because fluid elements evolve adiabatically, 
and the relationship between entropy and density depends upon 7 , the stability criterion 
in general depends upon the compressibility of the fluid. Notice that in the limit gL p <C 
c^, where L p = (9 In p/dr) -1 is the density gradient length scale, the growth rate is 
independent of 7 and reduces to the short-wavelength limit of the classical Rayleigh- 
Thylor growth rate when the effects of density gradient stabilization are included (e.g. 


Betti et all I1998T) . For the hot-spot configuration described above, it can be seen from 


Fig. [T] that the entire gas region and most of the shell is buoyantly unstable during 
deceleration. 

For a spherically symmetric mean flow (v = v r r ), perturbations can be decomposed 
into spherical harmonics Y) m and the resulting equations are 


dp' 

dt 


1_d_ 

r 2 dr 


1 d 


() P' ~ ) + 


1 ( 1 - 


1 ) , 
- PV_ L, 


(3.2) 


dv' r 

dt 

dv' ± 

dt 


dv r 




dr 


— v± 

r 


1 dp 
p 2 dr 


P 


P_ 

) 

pr 


ds' 

dt 


1 dp' 
p dr ’ 

ds . 

= “^ 


(3.3) 


(3.4) 


where a prime denotes a fluctuation, vf is the component of the velocity fluctuation 
parallel to ’V Y/ m , s' = p' / p — jp' /p is the entropy fluctuation, and i is a Legendre mode 
number ( Lai fc Goldrcichll2000ll . The third velocity component decouples from the other 


variables and will not be considered here. Short-wavelength, low-frequency fluctuations 
satisfy V • v' ss 0 and s' ss —7 p' / p (the Boussinesq approximation for buoyancy-driven 
flows; Johnson fc Gammid 2005 1. This reduces the perturbed continuity equation to the 
incompressive condition. It is important to note that the gas remains compressible; the 
incompressive condition applies only to the fluctuations. 

Notice from expression (13.11) that the growth rate of local modes in an equilibrium con¬ 
figuration is largest for modes with k r k±. This suggests the additional approximation 
of restricting the analysis to modes with d/dr -C £/r, i.e. modes that are sufficiently 
elongated in the radial direction that their radial variation can be neglected. I shall re¬ 
fer to this (along with the incompressive condition) as the oblate limit, and the modes 
thus isolated as oblate modes. This approximation considerably simplifies the analysis 
while retaining the essential physics. In the oblate limit, which amounts to neglecting the 
pressure perturbation in equation (13.31) , the linear equations reduce to the set of coupled 
ordinary differential equations 


dv' r 

dt 


dv r 




dr 


j_dp , 

p 2 dr P 


ds' 

dt 


ds 




dr 


(3.5) 


with s' = —7 p'/p . The first term on the right-hand side of the perturbed Euler equation 
represents angular momentum conservation: vortical modes spin faster (slower) under 
compression (expansion). The second term represents the baroclinic production of vor- 
ticity: entropy fluctuations spin up due to the baroclinic torque applied to them by the 
mean pressure gradient. The perturbed entropy equation represents the fact that, in 
the presence of a mean entropy gradient, entropy fluctuations evolve to compensate for 
mean entropy changes. An entropy fluctuation that moves up the mean entropy gradient, 
for example (to a region of higher mean entropy), must decrease in magnitude in order 
for total entropy to be conserved. This can increase or decrease the magnitude of the 
baroclinic term in the perturbed Euler equation. 

One can see the connection between the perturbed Euler equation in (13.51) and the 
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vorticity equation as follows. The dominant vorticity component for oblate modes is the 
one perpendicular to both r and Vyj m ; in axisymmetry, it is the only non-zero component 
and is given by 


/ v'j_ — v' r 9v' ± \ dPe(cosO) ~ v' r dPe(cosO) 
\ r dr ) 86 r 80 


(3.6) 


where the approximation is valid for oblate modes and Pg is a Legendre polynomial. The 
non-radial velocity of an oblate eddy is much smaller than its radial velocity fluctuation 
(this follows from V • v' « 0 and </r> d/dr), which implies that its vorticity is do mi¬ 
nated by radial motion in a frame moving with the mean flow. Appendix A of Johnson 
( 20141 ) derives the perturbed vorticity equation under the Boussinesq approximation; in 
spherical geometry and axisymmetry, the 4> component of expression (A4) of J ohnso n 
( 20141) (there is a sign error in front of the baroclinic term in that reference) is 


/ v r dv r \ 1 dpdPe(cosO) , 

dt \ r + dr ) ^ p 2 r dr dO ^ 


Using the approximation in (13.61) . one can readily show that equation (13.71) is equivalent 
to the perturbed Euler equation in (13.51) . 

Using the mean flow expression din (pr 2 ) /dt = — dv r /dr and recalling that the La- 
grangian derivative commutes with a derivatives but not r derivatives, equations (13.51) 
can be combined to give 


d 2 s' n dv r ds' , t9 , 
-d^ +2 ^ +Ns =°- 


(3.8) 


The only assumptions regarding the mean flow that have been made in deriving this 
equation are an ideal-gas equation of state and radial adiabatic flow. This equation 
therefore governs oblate modes in any flow satisfying these conditions. It is shown in 
appendix A that equation (13.81) also applies to flows with cylindrical and planar symmetry. 
For v r = 0, it reduces to the dispersion relation (13.11) with k r = 0. 

For Kidder implosions, the solutions to equation (13.81) have the form s' oc exp (— i j ui dt ), 
with w 2 = N 2 ; this demonstrates that the stability of Kidder implosions is governed by 
the Schwarzschild criterion. The full solution is given by 


where 


t 

P 


Cs 


= — cosh tp K -I- 1 r0 -- sinh cf) K , 

Po 

A cosh + 4 iWEM sinh , 

CsO Po CsO 





dt' 


o 1 ± ( t'/t c ) 


(3.9) 


(3.10) 


(3.11) 


and L s o = roi ( dso/da) 1 is the entropy gradient length scale at t = 0. In these expres¬ 
sions, all quantities with a zero subscript can be regarded as functions of a, since for 
oblate modes each fluid element evolves with time independently of all the others. For 
Nq = 0, entropy fluctuations are conserved and v' r = v' r0 /h (for p' Q = 0), reflecting the 
conservation of angular momentum as the gas is compressed. 

An implosion takes place over a finite time scale, such that the impact of unstable 
fluctuations on the mean flow depends upon both their initial amplitude and how much 
they are amplified over the course of the implosion. Expression (13.91) or (13.101) can be 
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Figure 2. Growth factor as a function of mode number for the profile shown in Fig. [I] and a 
stagnating Kidder implosion/explosion (with C = 10), showing results from the full linear code 
(solid) and the analytical oblate limit (dotted). 



t/t c 

Figure 3. Peak vorticity as a function of time for the profile given by expressions (1231 (1231 
(with <7 = 0.25 and a p = 1), a stagnating Kidder implosion/explosion (with C = 10), and 
i = 30. Shown are results from two-dimensional Kull calculations at several grid resolutions 
(solid, labeled by zones-per-wavelength), along with results from the linear code (dotted). 
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used to obtain an estimate of the growth factor when Nq < 0 : 

{p'/p)~ ■ ( torq l JV o*cl 


g = 


/final 


(P' / 7*) initial 


(2cy 

e 7r|Af 0 t c | 


(accelerating) 
(decelerating). 


(3.12) 


The growth factor for the decelerating solution takes into account the explosion phase 
that follows the implosion phase in the Kidder stagnating solution (i.e. h evolving from 
C to 1 and then back to C), since the gas is unstable both before and after stagnation. 
The same expressions apply to the amplification of vorticity fluctuations. 

Growth factors for a representative fluid element (a = 0.7) undergoing a decelerating 
Kidder implosion/explosion are shown in Fig. [5] as a function of Legendre mode number 
l. These results were obtained with a Lagrangian code that solves the full set of linear 
equations m m , where the radial profile of the initial entropy fluctuation was a 
Gaussian with a width 10% of the outer radius. It can be seen from Figure [2] that oblate 
modes are the fastest-growing modes and that fairly large values of £ are required to 
reach the oblate limit {£ > 100). Figure [3] shows the evolution of the peak vorticity 
for £ = 30 in both the line ar code a nd a two-dimensional version of the Lagrangian 
hydrodynamics code Kull 1 Rathkopj 2000h at various grid resolutions. These results 
indicate that a resolution of ^30 zones per wavelength is required to accurately capture 
the amplification of vorticity under an implosion in a numerical simulation. 

For implosions whose trajectory is given by expression (12.41) . the solutions to equation 
(|3.8I) again have the form s' oc exp (—ifeo dt ), where here 


2 , 5-3 7 


iDu> — N 2 = 0 


and D = dv r /dr is the dilatation. These solutions can equivalently be expressed as 
s' oc , where 


-m, 

where Ri c = N 2 / D 2 is a compressive Richardson number. It can be seen that instability 
(wi > 0, where w; is the imaginary part of u>) corresponds to /3 r < 0, where j3 T is the 
real part of the negative branch of /3 (for an implosion, a quantity that grows with time 
decays with h). 

For 7 = 5/3, /3 = ±v/— Ri c = ±\/—Nq t 2 , which is real for Nq < 0. For a monatomic 
gas, then, the stability of these implosions is also governed by the Schwarzschild criterion. 
For 7 < 5/3, however, the system is unconditionally unstable as a result of the additional 
growth due to compression. For A 2 < 0, i.e. Ri c > ([5 — 3y] /4) , the fluctuations have 
an oscillatory component with an amplitude that increases slowly with time. For A 2 > 0, 
the growth is purely a power law in time, and increases with decreasing Ri c . The critical 
Ri c above which oscillations appear varies between 0 for 7 = 5/3 and 1/4 for 7 = 1. This 
is reminiscent of the Richardson criterion for the stability of a stratified shear flow ( Milesl 


1961 : Chimonaslll970ll . with shear replaced by dilatation. In this case the transition is not 


between stability and instability, but rather simply between the presence and absence of 
oscillations. 

The complete linear solution for oblate modes under implosions satisfying expression 

fljgP is 


= h 


Po 

Po 


5 _ 37 

cosh 6 -\ -——^ sinh 6 

4A 


u r 0 


2Ri c 


Vro 3 [7 - 1]A 


sinh <f> 


(3.13) 
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V -^ = h- 


VrO 

V r Q 


cosh cj) — ° ^ sinh </> 


Po 3[7~ 1] 


Po 


2A 


sinh (j> 


(3.14) 


where <f> = VA 2 In h. For A 2 < 0, these solutions become oscillatory, and it is straight¬ 
forward to show that the amplitu de of the os cillations is consistent with the conservation 
of wave action in the fluid frame ( Whitham 1965ll . i.e. 


dt \ w r / 2 


1 .- 9^72 
r 5 


72 + 2^ 2 e 


where £ is the specific energy of the fluctuations (kinetic plus potential), = — L s s' is 
the fluctuating radial fluid displacement, w r is the real part of w, and an overbar denotes 
an angular average. An estimate of the growth factor for the fluctuations in this case can 
be made from (13.131) or (13.141) : 


C~ 


C<V-Ric 

C (5-3 7 )/2 

^(5-371/4 


for Ri c <C — 1 
for Ri c = 0 
for Ri c > 0. 


(3.15) 


4. Summary and discussion 

By isolating incompressive modes that are elongated in the direction of the mean flow 
(oblate modes), the following conclusions are drawn regarding the stability of homol¬ 
ogously imploding ideal gases: (i) monatomic gases are governed by the Schwarzschild 
criterion for buoyant stability; (ii) owing to the time-dependent nature of the mean flow, 
the growth is power law rather than exponential in time; (iii) additional growth occurs 
due to the conservation of angular momentum as vortices are compressed; and (iv) gases 
with 7 < 5/3 are weakly unstable due to this additional growth mechanism even when 
the f low is otherwise buoyantly stabl e. As point ed out bv ICao fc Loul ( 20091 ) . the reason 
that Goldreich h Webeil (1980) and Chu (119961) do not find instability is that they as¬ 
sume an isentropic background and ignore vorticity and entropy fluctuations; either of 
these assumptions precludes the development of buoyancy instability. 

The short-wavelength nature of the most unstable modes coupled with the compression 
of the mean flow makes this instability challenging to capture in numerical calculations. 
Fig. [3] indicates that accurately capturing the growth of vorticity under an implosion 
requires a resolution of the order of 30 zones per wavelength; this translates to an an¬ 
gular resolution of 360°/(30£) = \2° /£. Capturing the fastest-growing modes (£ > 100) 
therefore requires an angular resolution < 0.1°. For an Eulerian code, high convergence 
ratios impose severe constraints on the resolution in the direction of the mean flow. Seed¬ 
ing buoyancy instability in the gas can occur in two ways: vorticity fluctuations can be 
transported there by shocks rippled from drive asymmetry or interfacial perturbations, 
or small-amplitude entropy fluctuations can be present initially in the gas. Even a suf¬ 
ficiently resolved calculation that is initialized without ambient density or temperature 
fluctuations may not capture this instability properly. 

Finally, for ICF implosions and sonolumi nescent bubbles, the growth at large £ is likely 


to be reduced by conduction and viscosity (jAtzeni fe Mever-ter-Vehnll2004t IWeber et al 


12014 . These effects have been neglected here for several reasons, (i) As with classical 
stability analyses, a better understanding is gained if the effects driving the instability 
are isolated first, and stabilizing effects are added afterwards, (ii) There are uncertain¬ 
ties associated with conduction and viscosity models, and it is useful to know what 
instabilities are lurking in the background in their absence, (iii) The analytic results 
obtained here in the adiabatic limit provide physical insight and are useful for code ver- 
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ification. Self-similar solut i ons have been ob tained previously with non-adiabatic effects 
( Basko fc Murakami 1998 : Sanz et al. 2005), although the adiabatic profiles in Fig. [1] 
are remarkab ly s imila r to t he actu al p rofile s in an ICF implosion (compare Fig. |T] with 
Fig. 3.11 of Atzeni fc Mever-ter-Vehn 2004). Dissipation affects small scales; in reality 
the amplification in Fig. [2] will fall off at high mode numbers. At the same time, Fig. [2] 
indicates that potentially significant growth can occur even for moderate t. Definitive 
conclusions regarding the application of this analysis to ICF and sonoluminescent bub¬ 
bles will require a more faithful representation of both the mean flow and the dissipation. 
This, as well as application to astrophysical gases, will be pursued in future studies. 


I thank Dan Clark, Omar Hurricane, Karnig Mikaelian, Oleg Schilling, and the referees 
for their comments. This work was performed under the auspices of Lawrence Livermore 
National Security, LLC, (LLNS) under Contract number DE-AC52-07NA27344. 


Appendix A. Results for cylindrical and planar symmetry 


A complete set of homologous solutions can also be derived for flows with cylindrical 
and planar symmetry. In general, mass and entropy conservation for homologous adia¬ 
batic flow imply pit, a) = po(a)/h n and p{t,a) = po(a)/h nj , where n = 1, 2 or 3 for 
planar, cylindrical or spherical symmetry, and the time-dependent portion of ( 12.211 is 
h n ~ < -n+ 1 h = it,~ 2 . All of the above results for the mean spatial profiles remain valid if 
r is interpreted as the spatial coordinate in the direction of the mean flow. Expression 
(12.31) for the Kidder implosion remains valid for 7 = 1 + 2/n, and the generalization of 

m is 


h{t) 



n'y — n + 2 
2 


t 


n'y — nt c 


2 

717 — 71 + 2 


valid for hot c = — \j2Hny — n). Both the buoyancy frequency and the dilatation have 
the time dependence N oc D oc h/ n ~ ni ~ 2 ^ 2 . 

Equations (13.51) for the perturbations remain valid in all three geometries, and imply 


d 2 s' ds dv' r d f pr n 1 ds \ , 

dt 2 dr dt dt \porg ~ 1 dr 0 ) Vr ' 

where r$ = r(0, a). Using the mean flow relation din (pr 11-1 ) /dt = —dv r /dr, this leads 
to equation (13.81) . The solutions for a Kidder implosion, equations (13.91) and (13.101) . are 
therefore valid for all three geometries (provided 7 = 1 + 2/n). Expressions (13.131) and 
(13.14[) generalize to 



Po 

PQ 


, , 2 + n — n 7 . 

cosh 0 H---sinh < 

4A 


v'rp 2Ri c 
Vr0 n [7 - 1]A 



(Al) 



f v ' r0 

\VrO 


cosh cj) 


2 + n — n'y 
4A 


sinh (f> 


Po n[l - 1] 
po 2A 


sinh (j )) , 


(A 2) 


with A 2 = (2 + n —n7) 2 /16 —Ri c , so that oscillations appear for Ri c > (2 + n — n'y) 2 /16. 
For Ri c = 0, p' 0 = 0 and 7 = 5/3, expression (IA 21) reduces to v' r /v r = iv' r0 /v r o) h ra / 3-1 , 
which implies that the growth due to compression in a monatomic gas can impact the 
mean flow for planar and cylindrical implosions, but not for spherical implosions. In the 
latter case, the velocity fluctuations grow at the same rate as the mean and therefore 
never become large enough to drain energy from it. 
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